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STUDY QUESTION: Can high-resolution array comparative genomic hybridization (CGH) analysis of DNA samples from women with 
primary ovarian insufficiency (POI) improve the diagnosis of the condition and identify novel candidate genes for POI? 

SUMMARY ANSWER: A mutation affecting the regulatory region of growth differentiation factor 9 (GDF9) was identified for the first time 
together with several novel candidate genes for POI. 

WHAT IS KNOWN ALREADY: Most patients with POI do not receive a molecular diagnosis despite a significant genetic component in the 
pathogenesis. 

STUDY DESIGN, SIZE, DURATION: We performed a case-control study. Twenty-six patients were analyzed by array CGH for identi- 
fication of copy number variants. Novel changes were investigated in 95 controls and in a separate population of 28 additional patients with POI. 
The experimental procedures were performed during a I -year period. 

PARTICIPANTS/MATERIALS, SETTING, METHODS: DNA samples from 26 patients with POI were analyzed by a customized IM 
array-CGH platform with whole genome coverage and probe enrichment targeting 78 genes in sex development. By PCR amplification and 
sequencing, the breakpoint of an identified partial GDF9 gene duplication was characterized. A multiplex ligation-dependent probe amplification 
(MLPA) probe set for specific identification of deletions/duplications affecting GDF9 was developed. An MLPA probe set for the identification of 
additional cases or controls carrying novel candidate regions identified by array-CGH was developed. Sequencing of three candidate genes was 
performed. 

MAIN RESULTS AND THE ROLE OF CHANCE: Eleven unique copy number changes were identified in atotal of I I patients, including a 
tandem duplication of 475 bp, containing part of the GDF9gene promoter region. The duplicated region contains three NOBOX-binding elements 
and an E-box, important for GDF9 gene regulation. This aberration is likely causative of POI. Fifty-four patients were investigated for copy number 
changes within GDF9, but no additional cases were found. Ten aberrations constituting novel candidate regions were detected, including a second 
DNAH6 deletion in a patient with POI. Other identified candidate genes were TSPYL6, SMARCCI , CSPG5 and ZFR2. 
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LIMITATIONS, REASONS FOR CAUTION! This is a descriptive study and no functional experiments were performed. 

WIDER IMPLICATIONS OF THE FINDINGS: The study illustrates the importance of analyzing small copy number changes in addition to 
sequence alterations in the genetic investigation of patients with POL Also, promoter regions should be included in the investigation. 
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Introduction 

A clinical diagnosis of primary ovarian insufficiency (POI) is defined by the 
presence of primary or secondary amenorrhea (PA orSA) and hypergo- 
nadotropic hypogonadism, occurring before the age of 40 (Nelson, 
2009; De Vos et al., 2010). This condition has previously also been 
termed pre mature ovarian failure in patients with SA and gonadal dysgen- 
esis in PA. 

POI pathogenesis can be divided into follicle dysfunction and 
follicle depletion. Follicle depletion can be caused by a small initial 
germ cell count or a rapid germ cell loss. External factors, such as 
surgery, irradiation and infection as well as autoimmune and meta- 
bolic disorders also lead to POI. Most cases however remain 
idiopathic, and there is a significant genetic component (Simpson and 
Rajkovic, 1 999; Persani et al., 20 1 0). The most common are chromo- 
somal aberrations including Turner's syndrome (Baronchelli et al., 
20 1 2). Premutation of the FMRI gene is responsible for 2-5% of spor- 
adic cases with POI presenting with SA (Wittenberger e£ al., 2007; 
De Vos et a/., 2010). 

Other known genetic causes are mutations in genes for oocyte- 
secreted factors, GDF9 (growth differentiation factor 9) and BMP 1 5 
(Dixit eto/., 2006; Kovanci et al., 2007), for oocyte nuclear transcription 
factors NOBOX and FIGLA (Qin etai, 2007; Zhao eto/., 2008) and for the 
hormone receptors FSHR and NR5A I (Aittomaki eto/., 1 995; Janse eto/., 
20 1 2). However, the majority (90%) of patients with POI do not receive a 
molecular diagnosis (Nelson, 2009) and it is likely that many causative 
genes still remain unknown. 

By genome-wide detection of small genomic rearrangements, it is 
possible to identify novel genes involved in ovarian development, 
both on the X chromosome and on autosomes. The aim of this study 
was to identify novel candidate genes for ovarian development and 
function by investigating submicroscopic genetic imbalances in patients 
with POI. A I M array-CGH (comparative genomic hybridization) 
platform with 2.2 kb probe spacing for whole genome coverage and 
enrichment targeting 78 genes involved in gonadal development 
was used. At present, all studies reporting copy number variants 
(CN Vs) in patients with POI by array-CGH have used lower resolution 
platforms (Aboura et al., 2009; Ledig et al., 2010; Liao et al., 201 I). 
Studies using whole genome SNP array (McGuire et al., 201 I) and 
targeted analysis of the X chromosome using array-based techniques 
have also been reported (Dudding et al., 2010; Quilter et al., 2010; 
Knauffet al., 201 I). 



Materials and Methods 

Patients 

Fifty-four patients with POI (23 with PA, 3 I with SA) referred to the Clinical 
Genetic Laboratory of Karolinska University Hospital Stockholm, Sweden, 
were included in the study. Samples from 26 patients ( 1 7 with PA and 9 
with SA) were available for array-CGH experiments. 

The diagnosis of POI was set by PA or SA in a girl with female external geni- 
talia as well as internal Mullerian structures (uterus) and hypergonadotropic 
hypogonadism (FSH > 30 IU/1), determined twice. In the group with SA the 
overall median age at diagnosis was 22 years (12-37) (Table I). The nine 
patients with SA included in the array-CGH experiments were diagnosed 
at a median age of I 6 ( I 3-22) years. 

One patient with PA has several affected sisters with POI who did not want 
to be included in the study. Among patients with SA, four have mothers with 
POI, one a sister with POI and one an affected paternal aunt, who is included 
in the study. 

None of the patients had undergone previous ovarian surgery, chemother- 
apy or radiotherapy or presented with syndromic forms of POI. Five patients 
presented with hypothyroidism with anti-thyroid autoantibodies, after POI 
diagnosis in four with SA, and at I I years of age in one patient with PA. 
None had autoantibodies against adrenal cortex or 2 1 -hydroxylase protein. 

Nine patients with POI underwent laparoscopy with diagnostic ovarian bi- 
opsies. In four patients with PA, streak gonads with no visible follicles were 
found. Three patients with SA also had streak gonads with no visible follicles. 
One patient with SA at 23 had multiple primordial follicles and one patient 
with SA at I 5 had a few visible primordial and primary follicles. 



Table I Clinical overview of investigated patients. 



Clinical data 


PA 

(n = 23) 


SA 

(n = 3l) 


Median age at diagnosis, years (range) 


16(13-18) 


22(12-37) 


Median age at menarche, years (range) 




13 (10-15) 


Median FSH level at diagnosis, IU/1 


89 (39-150) 


80 (34- 1 55) 


(range) 3 






Swedish Caucasian, n (%) 


17(74) 


26 (84) 


Known heredity for POI, n (%) 


1 (4) 


6(19) 


Spontaneous conception before 




1 1 (35) 


diagnosis, n (%) 







PA, primary amenorrhea; SA, secondary amenorrhea. a Reference value: <30 IU/1. 
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Ethical approval 

All participants gave informed consent and the regional Ethics Committee at 
Karolinska Institutet, Sweden, approved the study with registration number 
2007/263-3 1 /2 and 20 1 I /276-32. 

Genetic investigation 

A 46, XX karyotype in peripheral blood was confirmed in all but in Patient 
24 in whom a balanced Robertsonian translocation identified (45, XX, der 
(13; 14)) was not considered causative. 

Sex chromosome mosaicism was excluded using fluorescent in situ hybrid- 
ization and DNA probes from chromosomes X and Y on peripheral blood 
smears, and when available, on touch preparations from gonadal tissue. 
Genetic investigations included sequencing of the NR5AI, BMP 1 5, FSHR, 
FIGLA, NOBOX and GDF9 genes. FMRI mutations were also excluded in 
patients with SA. 

Controls 

DNA samples from 95 healthy women were collected. All women were 
above the age of 40 and had given birth to at least one child. The exclusion 
criteria were previous egg donation, in vitro fertilization, infertility treatment 
or menopause before the age of 40. 

DNA extraction 

DNA was extracted from peripheral blood lymphocytes. Some samples 
were further purified using the QiAmp DNA minikit (QIAGEN, Sweden) 
to achieve acceptable quality values for array-CGH analysis. 

Array-CGH 

Twenty-six unrelated patients with POI were analyzed using a customized 
I M oligomarker array-CGH platform developed at Oxford Gene Technol- 
ogy (OGT, UK). In addition to whole genome coverage, the platform is 
enriched with probes targeting 78 genes implicated in sex development (Sup- 
plementary data, Table SI). The theoretical average probe resolution is 
2.2 kb. 

Preparation of labeled DNA and subsequent hybridization were per- 
formed according to the Agilent oligonucleotide array-based CGH for 
genomic DNA analysis' protocol (v6.2) and as previously described 
(Norling et al., 2013). A commercial DNA sample with pooled human 
genomic DNA from 1 0 female controls (Promega, Sweden) was used as ref- 
erence DNA. 

Data was analyzed using the Cytosure interpret software v3.4.3 (OGT). 
Circular binary segmentation analysis to detect copy number changes was 
performed using the following parameters: minimum probe count: 5, thresh- 
old for gains: 0.35, threshold for losses: 0.65, chromosome average method: 
median segment. 

By comparing data with the online database of genomic variation (DGV) 
version 10 released Nov 2010 (lafrate et al., 2004), common CNVs found 
in healthy control samples were excluded from further investigation. This 
exclusion was done with caution, as phenotypic effects of genetic imbalances 
affecting sexual development can be dependent on chromosomal sex. This 
information is not always available for controls included in the DGV. Aberra- 
tions only partially overlapping with rare reported CNVs were not excluded. 
Small intronic variations and aberrations not affecting genes were excluded 
after verification that they were not located just upstream or downstream 
a known gene causing disorders of sex development as a positional effect 
should then be considered. 



Multiplex ligation-dependent probe 
amplification 

Multiplex ligation-dependent probe amplification (MLPA) (Schouten et al., 
2002) probes were designed according to the recommendations by Stern 
et al. (2004) and combined in different probe sets (Supplementary data, 
Table Sit). The reference probes and the pilot probe used have been 
described earlier (Barbaro et al., 2008). 

MLPA reactions were performed according to the EKI reagent kit 
(MRC-Holland, the Netherlands) recommendations using 200 ng of DNA 
and the in-house designed probe set, as previously described (Barbaro 
et al., 2008). The reference DNA sample used for array-CGH was used as 
a control in all runs, along with at least another control DNA. 

Confirmation of array-CGH findings 

Copy number variations left after exclusion criteria were confirmed by MLPA 
usingtwo probes peraberration. When available, parental samples were ana- 
lyzed forcontrol of inheritance. Trace data were exported and analyzed in an 
Excel 2007 spreadsheet. Each sample's peak heights were normalized to the 
average peak height of the reference probes and subsequently normalized to 
the average peak height of the control samples. The analysis was accepted if 
the reference probe ratios were between 0.8 and 1 .2. Threshold values for 
deletions and duplications were set at 0.75 and 1 .25, respectively. 

MLPA for GDF9 

A specific synthetic probe set for MLPA analysis of the GDF9 gene was 
designed. Two probes were placed in each exon, two upstream GDF9 and 
one probe downstream (Fig. I E). The probe set was validated for consistency 
by analysis of 1 0 healthy controls, with a standard deviation <0. 1 for each 
probe. Trace data were analyzed using the GeneMarker v 1 .90 (Soft Genetics, 
USA) software using internal control probe normalization, quantification by 
peak height and ratio threshold values for deletion and duplication were 0.75 
and 1 .3, respectively. 

Patient cohort and control screening 

A specific synthetic probe set for MLPA analysis of array findings was 
designed, using one probe per each novel candidate copy number variation 
identified by array-CGH experiments. In addition, one probe within the 
keratin-associated protein (KRTAP)2-I and the KRTAP2-2 gene, respectively, 
were included. The probe set was validated for consistency. Trace data were 
analyzed with GeneMarker v 1 .90. 

Sequencing of TSPYL6, KRTAP2-3 and 
KRTAP2-4 

The single exon genes TSPYL6, KRTAP2-3 and KRTAP2-4, including part of the 
UTRs, were amplified by PCR using DyNAzyme EXT polymerase. Primers 
were designed using Primer3 software (Koressaar and Remm, 2007; Unter- 
gasseret al., 2012) (Supplementary data, Table SIN). 

PCR products were cleaned and sequenced using the ABI BigDye Termin- 
ator v3. 1 kit (Applied Biosystems, Sweden). Electropherograms were ana- 
lyzed against the reference sequences: NM_00 1 003937.2 for TSPYL6, 
NM_00l 165252.1 for KRTAP2-3 and NM_033I84.3 for KRTAP2-4, using 
the SeqScape v2.5 program (Applied Biosystems). 

Characterization of the GDF9 duplication 

Primers for the amplification and sequencing of the duplication junction were 
designed using Primer 3 software (Rozen and Skaletsky, 2000) (Supplemen- 
tary data, Table Sill) with the hypothesis of a tandem duplication. 
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In silico analysis 

For all identified candidate genes the following databases were searched for 
information. Relevant data are presented in the discussion section. 

NCBI (http://www.ncbi.nlm.nih.gov/) including PubMed, UCSC (http:// 
genome.ucsc.edu/) (Dreszer et al., 2012), GeneCards (http://www. 
genecards.org/), The Human Protein Atlas (http://www.proteinatlas.org/) 
(Uhlen et al., 201 0), Gene expression profiles during sex determination by 
Dr Serge Nef (http://nef.unige.ch/microarrays.php) (Nef et al., 2005), DE- 
CIPHER (Database of Chromosomal Imbalance and Phenotype in Humans 
Using Ensembl Resources) (http://decipher.sanger.ac.uk/), LifeMap Data- 
base of Embryonic development, Stem Cell Research and Regenerative 
Medicine (http://discovery.lifemapsc.com) (Edgar et al., 2013). 

Polyphen2 was used to evaluate amino acid substitutions (http://genetics. 
bwh.harvard.edu/ pph2). 

Results 

Array-CGH 

A total of 1 720 aberrations were detected by array-CGH analysis of 26 
samples with an average of 66 changes per patient (ranging from 29 to 
142). After exclusion criteria, 1 3 copy number changes remained, 6 in 
the PA group and 7 in the SA (Table II). Two pairs of patients carried 
the same aberration and one patient carried 3 changes, leaving I I 
unique copy number changes identified in a total of I I patients. One ab- 
erration affected the GDF9 gene, already associated with POI, the other 
1 0 changes unraveled novel candidate regions. All aberrations were con- 
firmed by MLPA. When possible, inheritance pattern was investigated. 

The entire patient cohort of 54 patients with POI and 95 healthy con- 
trols were screened by MLPA for additional cases with copy number 
changes within the novel candidate regions. One patient and two con- 
trols were heterozygous for KRTAP2-3 and KRTAP2-4 deletion. One 
control was hemizygous for the KRTAP2-2 gene. No other patient or 
control was found to carry any other identified aberration. 

GDF9 duplication investigation and GDF9 
MLPA analysis 

A duplication between 475 and 1 729 bp within the GDF9 gene was iden- 
tified by array-CGH in Patient 20. We amplified the duplication junction 
(Fig. I D) and by sequencing, we confirmed a tandem head-to-tail dupli- 
cation of 475 bp containing the first part of the 5' UTR and a short up- 
stream DNA sequence (Chr5.hg 1 9:g. 1 32,20 1 , 1 70_ 1 32,20 1 ,644dup). 

Analysis of the entire patient cohort by an MLPA probe set for the 
identification of deletions/duplications within or encompassing GDF9 
did not identify any additional case. 



Sequencing 

The KRTAP2-3 and KRTAP2-4 genes were sequenced in the two hemizy- 
gous patients. No mutations were identified in KRTAP2-4, while the vari- 
ation c.2I8G>A (p.Cys73Tyr), that is reported in dbSNP as 
rsl 13397060, was identified. KRTAP2-3 sequencing was extended to 
parental samples and the entire patient-control cohorts. One more 
patient and one control were homozygous (A/A) for this substitution. 

Sequencing of the TSPYL6 gene in the entire patient cohort revealed no 
mutations. 

Discussion 

Out of a group of 26 patients with POI, we have identified I I unique 
novel aberrations in I I patients, using a customized high-resolution 
array-CGH platform. Our detection rate is comparable to Ledig et al. 
(2010), and higher than studies using platforms with lower resolution 
(Aboura et al., 2009; Liao et al., 20 1 I ; McGuire et al., 20 1 I ). 

One change affects GDF9 and the other 1 0 involve new regions. Ma- 
ternally inherited changes, even though less likely pathogenic if assuming 
a haploinsufficiency mechanism, have not been excluded by default as 
autosomal recessive inheritance, and polygenic or environmental 
factor affecting phenotype penetrance cannot be excluded. 

GDF9 on 5q3 I . I 

A duplication affecting GDF9 was detected in Patient 20 with SA at 1 5 years 
of age (Fig. I ). This oocyte-secreted factor, expressed from the primary fol- 
licle stage, is necessary for normal folliculogenesis and fertility. Mutations 
affecting the open reading frame of the GDF9 gene have been associated 
with POI at a low frequency (Dixit et al., 2005; Laissue et al., 2006; Zhao 
eta/., 2007). No such mutations were identified in our patient cohort. 

Within the duplicated region there are three NOBOX-binding ele- 
ments (NBEs) and an E-box sequence. In mouse NBEs regulate Gdf9 ex- 
pression (Choi and Rajkovic, 2006), and the E-box sequence is a target 
for transcription factors necessary for the ovarian specific expression 
of Gdf9 (Yan et al., 2006). Both the NBEs and the E-box are conserved 
between mice and humans, as well as several other species. Thus, they 
are likely important also for the regulation of human GDF9 expression 
(Choi and Rajkovic, 2006; Yan etai, 2006). It is likely that the duplication 
causes altered GDF9 expression in the developing ovary, leadingto POI. 
This represents the first mutation identified affecting the regulatory 
region of GDF9. 

An enrichment of CNVs has been identified in promoter regions, sug- 
gesting that promoter sequences can be mildly unstable. This is a 



Figure I GDF9 duplication. (A) Representation of an array-CGH result from the Cytosure software. The duplicated segment is indicated by blue back- 
ground with a positive baseline offset. Blue dispersed dots represent oligomarkers. Blue arrows on bottom indicate gene location and purple smaller arrows 
represent exons. Black horizontal line on bottom represents location of duplication. Chromosome ideogram shown on top with the current segment mag- 
nification indicated by blue background color. (B) Schematic representation of the GDF9 gene. The ATG indicates the initiation of translation. Red dots 
represent N BE. Violet dot represents the E-box. Yellow arrow indicates the duplicated segment. Blackand blue arrows represent PCR primers for wild-type 
allele (wt) and for duplication breakpoint amplification, respectively. (C) Schematic representation of the duplication. The green line represents the PCR 
product containing the duplication junction. (D) Breakpoint characterization. Gel image of amplified PCR breakpoint fragment. P, patient; C, control; B, 
blank; L, ladder. Electropherograms with sequence of the breakpoint (in the middle), partially aligned with the wt allele (upper and lower electrophero- 
grams). (E) Representation from the UCSC genome browser, GRCh37/hg 1 9 assembly. Horizontal yellow arrow line indicates the duplicated segment. 
Locations of NBE (TAATTA, TAATTG) and E-box (CAGCTG) are shown as vertical red and purple lines, respectively. Vertical black boxes represent 
MLPA probes. 
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Table II Results from 26 patients, copy numbers shown after application of exclusion criteria. 



Pat. Age at Chr. Start End Size Probe Del/ Gene(s) Inheritance Other information 

no. amen. band (kb) count dup 



1 


PA 


2p 1 6.2 


54296546 


54398 1 00 


102 


39 


Del 


ACYP2, TSPYL6 


Not maternal 




2 


PA 


















Ovarian biopsy 


3 


PA 


I5q26.2 


92580470 


92619918 


39 


17 


Del 


MCTP2 


Maternal 




4 


PA 


3p2l.3l 


47588558 


47723275 


[ 35 


54 


Dup 


SMARCCI, 
CSPG5 


Maternal 


Mother breast cancer 4 1 
years 


5 
6 
7 


PA 
PA 
PA 
PA 

PA 
rA 

PA 
rA 


I5q26.2 


92580470 


92619918 


39 


17 


Del 


MCTP2 


Paternal 




8 

Q 

7 

1 U 


1 Zp 1 j.j j 


ZODOO/ 1 


LI 1 o/LO 


jZ 


1 D 


LJUp 




FVl ^ 1 

1 laLernai 


- 


1 1 
12 
13 
14 


PA 
PA 
PA 
PA 


I5q25.2 


81917995 


81997930 


8 


25 


Dup 


SH3GL3 


Maternal 


- 


15 
16 
17 


PA 
PA 
PA 


- 
— 


- 
— 


— 
— 


— 


- 
— 


— 


— 
— 


— 
— 


— 
— 


18 


1 6 years 




















19 


1 4 years 


1 9p 1 3.3 

Xp22.33 
I7q2l.2 


3688186 

1529221 
36465290 


38161 17 

1572694 
36473792 


■ 

9 


51 

17 
5 


Del 


RAX2, JJP3, 
ZFR2, MRPL54, 
APBA3, MATK 
P2RY8, ASM 71 
KRTAP2-3, 
KRTAP7-4 

t\t\ t f\t I- i 


Maternal 

Not maternal 
Not maternal 


Hypothyroidism ovarian 
biopsy 

- 


on 


1 5 years 


^■J 1 1 

jq J 1 . 1 


1 JZZZOZOO 


i imsz7dr\ 

1 j/Z/O/'rU 


U.D 


7 
1 


Pli ir% 

uup 


(jury 


M A 
IMA 




21 


22 years 


2pl 1.2 


84620455 


84791661 


171 


57 


Del 


DNAH6 


NA 


Paternal aunt POI 34 


22 


20 years 


I7q2l.2 


36465290 


36473792 


9 


5 


Del 


KRTAP2-3, 
KRTAP2-4* 


Mother 
hemizygous 


Hypothyroidism 
consanguinity mother 
with bilateral breast 
cancer 


23 


1 7 years 


2q34 


210703481 


210739467 


36 


14 


Del 


KANSLIL 


Maternal 




24 


21 years 


















45,XX,der(l3;l4) 
Ovarian biopsy 


25 


1 5 years 




















26 


1 3 years 


















Hypothyroidism with 
anti-TPO ovarian biopsy 



All coordinates are given using the NCBI36/hgl 8 build. Pat no., patient number; Age at amen., Age at amenorrhea; Chr. Band., chromosomal band; Del., deletion; Dup., duplication; PA, 
primary amenorrhea; NA, no available sample, 
homozygous deletion. 



mutations of TSPYLI (TSPY-Like I) have been found to cause dysgenesis 
of the testis and sudden infant death syndrome (MIM 608800). We 
considered TSPYL6 an interesting candidate gene for POI and 
sequenced the entire patient cohort, but no inactivating mutations 
were identified. 

ACYP2 encodes a muscle type form of acylphosphatase 2 and we con- 
sider it a less likely candidate for POI. The deletion occurs just upstream 
of C20RF73, an uncharacterized gene expressed in both ovaries and 
testis. A positional effect of the deletion on this gene expression 
cannot be excluded as a potential cause of POI. 



suggested evolutionary mechanism for changes in gene regulation 
(Conrad et a/., 2010). 

Candidate genes 

TSPY-Like 6, acylphosphatase 2, muscle type I and C20RF73 on 
2pl6.2 

The deletion, affecting TSPY-Like 6 (TSPYL6) and acylphosphatase 2, 
muscle type I (ACYP2), in Patient I was not inherited from the 
mother. The TSPYL6 gene function is unknown, but homozygous 
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Multiple C2 domains, transmembrane 2 on I5q26.2 
In two unrelated patients with PA (Patients 3 and 5), a 39 kb deletion was 
detected, removing a large part of intron I in one of the several isoforms 
of the multiple C2 domains, transmembrane 2 (MCTP2) gene 
(ENST00000543482). MCTP2 is a transmembrane protein that binds 
calcium in the absence of phospholipids. The deletion could affect spli- 
cing or regulation or have no effect at all. It is difficult to find an immediate 
connection between MCTP2 and POI with the available data. 

SWI /SNF-related matrix-associated, actin-dependent regulator of 
chromatin subfamily c member I and chondroitin sulfate proteoglycan 
5 on 3p2l.3l 

A 1 35 kb duplication on chromosome 3p2 1 .3 I was found in Patient 4 
with PA, spanning the SWI/SNF-related matrix-associated, actin- 
dependent regulator of chromatin subfamily c member I (SMARCCI) 
gene and at least the first two exons of the CSPG5 gene. SMARCC I is 
a member of the SWI/SNF family of proteins that regulate transcription 
by chromatin remodeling. It is highly expressed in both follicle and ovarian 
stroma cells. Chondroitin sulfate proteoglycan 5 (CSPG5) is a proteogly- 
can that may function as a differentiation factor. In mice, there is a sex dif- 
ferential expression profile of Cspg5 in Sfl+ somatic cells in developing 
gonads at embryonic day 12.5 (El 2.5) and El 3.5 (Nef etal., 2005) indi- 
cating a potential role in gonad differentiation. We consider both 
SMARCCI and CSPG5 possible candidate genes for POI. 

Calcium channel, voltage-dependent I type alpha I C subunit on 
1 2p 1 3.33 

Patient 9, with PA, has a maternally inherited duplication involving the 
calcium channel, voltage-dependent I type alpha I C subunit (CACNA I Q 
gene, likely not affecting expression as loss of function mutations causes 
autosomal dominant cardiac conduction abnormalities (Brugada syn- 
drome, MIM 61 1875; Timothy syndrome, MIM 601005). We do not 
consider CACNAICz candidate gene for POI. 

SH3-domain GRB2-like 3 on I5q25.2 

Patient 1 2, presenting with PA, carries a maternally inherited duplication 
of almost the entire sequence of the first intron of the SH3GL3 gene. 
SH3-domain GRB2-like 3 (SH3GL3) expression is abundant in testis 
(Ringstad et ai, 1 997), where it is believed important for spermatogen- 
esis (Li etal., 2009). Moderate expression in follicle cells is reported in the 
Human Protein Atlas. SH3GL3 interacts with the similar protein SH3GL2 
(SH3-domain GRB2-like-2) (Franceschini etal., 20 1 3) that in mice Sh3GI2 
has a sex differential expression profile at El 2.5 and El 3.5 (Nef et ai, 
2005). An association between SH3GL3 and POI should be further 
investigated. 

Dynein axonemal heavy chain 6 on 2pl 1.2 

In Patient 21 we have identified a 171 kb deletion that removes at 
least 52 of the 77 exons of the dynein axonemal heavy chain 6 
(DNAH6) gene encoding a heavy axonemal dynein chain. Dyneins are 
microtubule-associated motor protein complexes, important for ciliar 
and flagellar motility. There are heavy, light and intermediate chains. 
The mouse intermediate chain Dnaic2 (dynein, axonemal intermediate 
chain 2) has been found on the surface of oocytes in secondary and 
antral stages, but not in primordial or primary follicles suggesting a role 
for dyneins in ovarian development (Yang and Wu, 2008). 



Interestingly, in a previous array-CGH study a smaller but partially 
overlapping deletion of DNAH6 was detected in a patient with secondary 
amenorrhea (Fig. 2) (Ledig et ai, 20 1 0) strengthening the possible role 
for DNAH6 in POI. Another study has also identified the related 
DNAH5 as a candidate gene for POI (Aboura et ai, 2009). Mutations 
in DNAH5 are otherwise known to cause primary ciliary dyskinesia 
(PCD, MIM244400). We believe that DNAH6 is a very interesting candi- 
date gene for POI. 

KRTAP2-3 and KRTAP2-4 on 17 q2 1. 2 

A 9 kb deletion encompassing the KRTAP2-3 and KRTAP2-4 genes was 
identified in two unrelated patients with SA and by MLPA screening 
one additional hemizygous patient was found. 

We hypothesized an autosomal recessive causative mechanism, and 
by sequencing identified a potentially damaging amino acid substitution 
in the KRTAP2-3 gene in the hemizygous patients and in homozygous 
form in one more patient with POI. However, also one control was 
found to be homozygous for the substitution, contradicting KRTAP2-3 
mutations as causative for POI. It is possible that the identified 
KRTAP2-3 and KRTAP2 -4 deletion is a normal variant in the Swedish popu- 
lation as also two controls are hemizygous. 

Purinergic receptor P2Y G-protein coupled 8 and acetylserotonin 
O-methyltransferase-like on Xp22.33 

In addition to a deletion of KRTAP2-3 and KRTAP2-4, Patient 22 also 
carries a duplication in the pseudoautosomal region I on the X chromo- 
some. The duplication includes the last exon of purinergic receptor P2Y 
G-protein coupled 8 (P2RY8) and the first exon of the ASMTL gene. 
P2RY8 is a G-protein coupled receptor expressed in lymphocytes. Acet- 
ylserotonin O-methyltransferase-like (ASMTL) is a generally expressed 
enzyme, with high levels in ovarian follicle cells. We consider both 
P2RY8 and ASTML less likely candidate genes for POI. 

Del I9pl3.3 

A 1 28 kb deletion on chromosome 1 9 was also identified in Patient 22 
affecting six genes: TJP3, APBA3, MRPL54, RAX2, MATK and zinc finger 
RNA-binding protein 2 (ZFR2). 

Knock-out mouse models for TJP3 (Xu et at. , 2008), APBA3 (Hara et al. , 
20 1 I) and MATK (Lee etal., 2006) have been described without particu- 
lar phenotypes and with normal fertility excluding these genes as candi- 
dates for POI. 

MRPL54 encodes for a mitochondrial ribosome subunit necessary for 
the synthesis of mitochondrial proteins (Koc etal., 200 1 ). Defects of mito- 
chondrial protein synthesis lead to a more wide and severe phenotype than 
isolated POI. RAX2 heterozygous missense mutations cause age-related 
macular degeneration (MIM 610362) and cone-rod dystrophy type I I 
(MIM 6 1 038 1 ). We exclude MRPL54 and PAX2 as candidate genes for POI. 

ZFR2 has no known function, but is expressed in many adult cell types, 
including ovarian follicle and stroma cells, as well as in female gameto- 
cytes and germ cells (Edgar et ai, 2013). Ovarian biopsies from the 
patient show bilateral complete lack of oocytes and follicles and we con- 
sider ZFR2 an interesting candidate gene for POI. 

KAT8 regulatory NSL complex subunit I -like on 2q34 
In Patient 23, with SA at 1 7, we detected a 36 kb deletion removing the 
first coding exon of the KAT8 regulatory NSL complex subunit I -like 
(KANSLIL) gene on chromosome 2q34. KANSLIL, together with 
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Figure 1 DNAH6 deletion. (A) Representation of an array-CGH result from the Cytosure software. The duplicated segment is indicated by blue back- 
ground with a negative baseline offset. Blue dispersed dots represent oligomarkers. Blue arrows on bottom indicate gene location and purple vertical lines 
represent exons. The DNAH6 gene is shown with previous name DNHLI . The black horizontal line on bottom represents the location of deletion. Chromo- 
some ideogram shown on top with the current segment magnification indicated by blue background color. (B) Representation from the UCSC genome 
browser, GRCh37/hgl9 assembly. The green line represents the deletion detected in Patient 21. The orange line shows the deletion reported by 
Ledig et al. Vertical black lines represent MLPA probes. 



KANSLI on chromosome 1 7, are the human homologs of the Drosophila 
Waharan gene (Lone et al., 20 1 0). KANSL I is a regulator of KAT8 influ- 
encing gene expression through histone H4 lysine 1 6 acetylation 
(Mendjan eta/., 2006). Mutations in this gene cause 1 7q2 1 .3 1 microdele- 
tion syndrome (MIM 6I2452). No specific information on KANSLI L 
function is reported, but as an altered chromatin acetylation could pos- 
sibly affect developing germ cells and follicle formation, we consider 
KANSLI L a potential candidate gene for POL 

Conclusion 

Using array-CGH, we have identified the first mutation affecting the regu- 
latory region of GDF9 in a patient with SA, most likely pathogenic. Import- 
antly, the small, 475 bp duplication would not have been detected 
without the customized probe enrichment, making clear that platform 
resolution and probe targeting are crucial. We have also identified a 
second DNAH6 deletion in a patient with SA, corroborating DNAH6 as 
a candidate gene for POL In addition, we describe several novel candi- 
date genes such as TSPYL6, SMARCCI, CSPG5 and ZFR2. 



We consider array-CGH a powerful tool to identify candidate genes 
for POI and we recommend that patients without a molecular diagnosis 
should undergo analysis for copy number alterations. Soon next gener- 
ation sequencing will be available and able to provide reliable data on 
copy numbers, once this occurs whole genome sequencing should be 
considered for diagnostics of POI patients, as most known causative 
genes can either harbor small pathogenic sequence changes or be 
present in abnormal copy numbers. 

Supplementary data 

Supplementary data are available at http://humrep.oxfordjournals.org/. 
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